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ABSTRACT 

We use observations of cosmic microwave background anisotropics, supernova lu- 
minosities and the baryon acoustic oscillation signal in the galaxy distribution to 
constrain the cosmological parameters in a simple interacting dark energy model with 
a time- varying equation of state. Using a Monte Carlo Markov Chain technique we de- 
termine the posterior likelihoods. Constraints from the individual data sets are weak, 
but the combination of the three data sets confines the interaction constant Y to be less 
than 23% of the expansion rate of the Universe Hq] at 95% CL —0.23 < T / Hq < +0.15. 
The CMB acoustic peaks can be well fitted even if the interaction rate is much larger, 
but this requires a larger or smaller (depending on the sign of interaction) matter 
density today than in the non-interacting model. Due to this degeneracy between the 
matter density and the interaction rate, the only observable effect on the CMB is a 
larger or smaller integrated Sachs- Wolfe effect. While SN or BAO data alone do not 
set any direct constraints on the interaction, they exclude the models with very large 
matter density, and hence indirectly constrain the interaction rate when jointly anal- 
ysed with the CMB data. To enable the analysis described in this paper, we present in 
a companion paper arXiv:0907.4981 a new systematic analysis of the early radiation 
era solution to find the adiabatic initial conditions for the Boltzmann integration. 

Key words: cosmology:theory, cosmology:observations, cosmic microwave back- 
ground, cosmological parameters, dark matter, large-scale structure of Universe 



1 INTERACTING DARK ENERGY 

Dark energy and dark matter are the dominant sources in 
the 'standard' model for the evolution of the universe. Both 
are currently only detected via their gravitational effects, 
with an inevitable degeneracy between them (one requires 
a model to separate dark energy from dark matter). There 
could thus be an interaction between them that is consis- 
tent with current observational constraints. A dark sector 
interaction could also alleviate the 'coincidence' problem 
(why are the energy densities of the two components of the 
same order of magnitude today?). Furthermore, interacting 
dark energy exerts a no n- gravitational 'drag' on dark mat- 
ter, and thus can introduce new features to structure for- 



mation, including possibly a new large-scale bias ( |Amen- 
dola & Tocchini-Valentini 2002| and a violation by dark 
matter of the weak equivalence principle on cosmological 
scales ( |Bertolami et al.|2007 Koyama et al. 2009). 

The energy balance equations in the background are 



p c = —3Hp c + aQc , 

Pde = -3H(l + W d e)pde+aQde 



(1) 

Qde = ~Qc , (2) 



where a is the scale factor of the Universe, Ti = a /a is 
the conformal Hubble parameter, Wde = Pde/ pde is the dark 
energy equation of state parameter, a prime indicates deriva- 
tive with respect to conformal time r, and Q c is the rate of 
transfer to the dark matter density due to the interaction. 

Various forms for Q c have been put forward (see, 



e.g. |Wetterich| fl!995|); |Amendola| fl!999|); |Billyard fe Coley 



fl2000|); |Zimdahi fe Pavonl (| 20011 ; |Farrar fe Peebles| (|2004[ 



Chimento et al. (2003); Olivares et al. (2005); Koivisto 



fl2005|);|Sadjadi fe Alimohammadi| ([20061 ;|Guo et ai.|fl2007|); 



Chongchitnan 



et al. (2009) 



2008); iPereira & Jesus (2009k; IQuercellini et al 



Boehmer et al.| (|2QQ8|); |He fe Wang| (|2QQ8|); |Quartin et at 



Valiviita et al | ( 2008 ); |He et al.| (|2009|; |Bean et al 



(2009) ; |Corasaniti| ( 2008|); |Caldera-Cabral 



Gavela et al. (2009) 



(2008); 



(2008); 



Jackson et al. (2009)) 



All of these models are phenomenological. Some of them are 
constructed specifically for mathematical simplicity - for ex- 
ample, models in which Q c oc Tip. Rather than design the 
interaction to achieve a specific outcome, we prefer to start 
with a simple physical model, and then develop its predic- 
tions. 
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Figure 1. Marginalized likelihoods for the interacting model with WMAP, WMAP&ACBAR, and WMAP&SN&BAO data. The 
curves show 68% and 95% CL regions: dotted magenta/grey for WMAP, solid cyan/grey for WMAP&ACBAR, and solid black for 
WMAP&SN&BAO data. 



We consider models which are similar to simple mod- 



els of reheating (Turner 1983), of dark matter decay into 
radiation (Cen 2001), and of curvaton decay ( |Malik et al. 
2003) - i.e., where the interaction has the form of a decay 



of one species into another, with constant decay rate. Such 



et al. (2008) 



a model was introduced in Boehmer et al. (2008); Valiviita 



Qc 



-Tp c 



(3) 



where Y is the constant positive decay rate of dark matter 
to dark energy. Here, as well as in Valiviita et al. (2008), we 
include also the possibility Y < 0, i.e., allow also for energy 
transfer is from dark energy to dark matter. 



In Valiviita, Majerotto & Maartens (2008) we consid- 



ered the case of fluid dark energy with a constant equation 
of state parameter — 1 < Wde ^ —4/5, and found that a se- 
rious large-scale non-adiabatic instability affects this model 
in the early radiation dominated epoch. This instability is 
stronger the closer Wde is to —1. Phantom models, Wde < — 1, 
do not suffer from this instability, but we consider them to 
be unphysical. 

The instability is determined by the early-time value of 
Wde- for constant Wde there is no flexibility, but a variable 
Wde (as in quintessence models, see e.g. Lee et al. (2006)) 



should allow us to avoid the instability. We show that the 
models are viable if Wde > —4/5 at early times, while at late 
times, Wde ~ — 1. We represent Wde via the parametrization 
Wde = wq + w a (l — a) ( |Chevallier &; Pola rski 2001 Linder 



2003 ) , which we rewrite as 



w d( 



— woa + w e (l 



(4) 



where w e = wo + w a is the early-time value of Wde , while wo 
is the late-time value. This parametrization was originally 
developed to phenomenologically account for the possible 
time evolution of Wde up to redshifts of a few. At higher 
redshifts it describes the simplest possible model where Wde 
is a constant, namely w e . In particular, this parametrization 
works well for some classes of quintessence models Caldwell 
& Linder ( 2005 ) and it has the advantage of having a finite 
value at high redshifts. Using a parametrisation for Wde has 



the drawback that it cannot reproduce all possible models 
for a large redshift range. Furthermore, it is not possible to 
compute the speed of sound but it is necessary to assume 
a value for it. It might be interesting to look at different 
parametrisations of Wde or a definite scalar field model. This 
is left for future work. In this paper we demonstrate with 
the simple parametrization, Eq. Q, that a (suitably) time- 
varying Wde cures the interacting model, Eq. (|3|, from the 
early-time large-scale instability, and thus provides a viable 
cosmology. 

We perform a full Monte Carlo Markov Chain (MCMC) 
likelihood scan for the spatially flat interacting and non- 
interacting models, using a modification of the CAMB 
\ Lewis et al. 2000| Boltzmann code, to determine the 
best-fitting values of Y/Ho and the other cosmological pa- 
rameters, against Wilkinson Microwave Anisotropy Probe 



(WMAP) 5-years data (Komatsu et al. 2009), Supernovae 
la (SN) Union sample data (Kowalski et al. 2008k, and 



data points for the ratio of the sound horizon to a dis- 
tance measure at two different redshifts from baryon acous- 



tic oscillation (BAO) observations (Percival et al. 2007). 
The best-fitting models with various combinations of data 
are shown in Table |A2| on page [14] in Appendix [B] while 
Fig. [I] above summarizes our main findings for the pos- 
terior likelihoods. Our most stringent results for the in- 
teracting model result from the combined analysis of 
WMAP&SN&BAO, giving the following minimal 95% in- 
tervals: u h G (0.0212, 0.0241), u c G (0.859, 0.125), H Q G 
(63, 70), r G (0.057, 0.133), Y/H G (-0.23, +0.15), w e G 
(-0.80, -0.19), w G (-1.00, -0.63), n s G (0.937, 1.002), 
ln(10 10 A|) G (2.95,3.14), Q d eo G (0.648,0.767), AgeG 
(13.6, 14.3) Gyr. Description and prior ranges of these pa- 
rameters are given in Appendix |A2[ There, in Table Al, we 
also list the definitions of some of the symbols used in this 
paper. 

The key features of the constraints on the interacting 
model from data may be summarized as follows. 



• Any of the data sets alone (CMB, or SN, or BAO) would 
allow for large interaction: |T| could be even larger than 
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Table 1. The evolution of perturbations on super-Hubble scales with various values of the dark energy equation of state parameter in 
the radiation and matter dominated eras. "Adiabatic" means that it is possible to specify adiabatic initial conditions so that the total 
gauge invariant curvature perturbation £ stays constant on super-Hubble scales, and the evolution of all non-dark energy perturbations is 
the same as in the non-interacting case while the de perturbations behave differently. "Adiabatic (standard)" means that the behaviour 
of all perturbations at early times on super-Hubble scales is the same as in the non-interacting model. 



in the RD or MD era 


Radiation dominated era 


(RD) 


Matter dominated era (MD) 


Viable? 


w de < -1 


adiabatic 




adiabatic 


viable, but phantom 


-1 < w de < -4/5 


"blow-up" isocurvature g 


rowth 


"blow-up" isocurvature growth 


non-viable 


-4/5 SC w de < -2/3 


adiabatic 




isocurvature growth 


viable, if T small enough 


-2/3 SC w de < -1/2 


adiabatic (standard) 




adiabatic 


viable 


-1/2 <:w de < +1/3 


adiabatic (standard) 




adiabatic (standard) 


viable 



today's Hubble rate Ho. In the CMB the only hint from a 
large interaction rate is a modified integrated Sachs-Wolfe 
(ISW) effect. However, due to the cosmic variance the \ 2 is 
only mildly affected. 

• A large negative T fits the CMB TT and TE spectra 
equally well as the F = model, but a good fit requires 
a larger physical cold dark matter density today, cj c , (and 
hence a smaller Qdeo) and a smaller H . Negative V sup- 
presses the late ISW effect and hence the CMB fit is slightly 
better than in the V — case. Combining the CMB data 
with either SN or BAO or SN&BAO data, this improve- 
ment is cancelled by a worse fit to SN and BAO due to too 
little acceleration at low redshifts. 

• Models with large positive V fit the high-/ CMB TT and 
all TE data equally well as the r = model, but a good fit 
requires a smaller uj c (and hence a larger Qdeo) and a larger 
Ho. The SN and BAO data can be fitted better than in the 
r = case due to increased acceleration at small redshifts. 
However, this improvement is cancelled by a worse fit to the 
low-/ CMB TT spectrum due to a large late ISW effect. 

There are two critical features of the analysis of inter- 
acting models, which are not always properly accounted for 
in the literature: 

• The background energy transfer rate Q c does not in it- 
self determine the interaction in the perturbed universe: one 
must also specify the momentum transfer rate. We do this 
via a physical assumption, i.e., that the momentum transfer 
vanishes in the dark matter rest-frame, so that the energy- 
momentum transfer rate is given covariantly ( | Valiviita et al.| 
[2008] ) by 



-r> c (i + £ c 



(5) 



where u% is the dark matter 4- velocity, and 5 C = Sp c /p c is 
the cold dark matter (CDM) density contrast. 

• Adiabatic initial conditions in the presence of a dark 
sector interaction require a careful analysis of the early- 
radiation solution. We derive these initial conditions in the 



companion paper (Majerotto, Valiviita & Maartens 2009) 



by generalizing the methods of Doran et al. (2003) to the 
interacting case, extending our previous results ( |Valiviita| 
|et al.|[2QQ8| ) . The key results for the initial conditions and 
early-time perturbation evolution are reproduced in Table [I] 

We give here the first analysis of the CMB spectra 
and the first MCMC likelihood analysis for the interact- 
ing model |5]), using the perturbation equations and ini- 
tial conditions given in the companion paper (Majerotto, 



Valiviita & Maartens 2009). Cosmological perturbations of 
other interacting models have been investigated in |Amen-| 



dola et al. 



(2003); 



Koivisto 



Mainini k Bonometto (2007) 



( |2005|); IQlivares et al. (2006) 
'Bean et al.| (|2008|); (Verganil 



LTJ20I 



et al.|(|2009|;|Pettorino fe Baccigalupi (2008); Schafer|(|2008|); 
Schaefer et al.|(|2008|);|La Vacca fe Colombo|p008^;|He et al 
(1 2009 1); [Bean et al. (2008); Corasaniti| fl2008j) ; |Chongchitnan 
|2009|; I Jackson et al.|(|2009|);|Gavela et al.|(2009|);|La Vacca 



et al. (2009); He et al. (2009); Caldera-Cabral et al. (2009) 



(2009) 



He et al. (2009); Koyama et al. (2009); Kristiansen et al. 



2 PHENOMENOLOGY 

We have performed 7 MCMC runs for the spatially flat 
(Q = 1) interacting model with various data sets (WMAP, 
WMAP&ACBAR, SN, BAO, WMAP&SN, WMAP&BAO, 
and WMAP&SN&BAO). Here WMAP refers to t he 5-year 
temperature and polarization anisotropy data 



_Komatsu| 

et al. |2009] ), ACBAR to t he Arcminute Cosmolo gy Bolome- 
ter Array Receiver data ( |Reichardt et al.||2009|), SN to the 
Union Supernovae la sa mple ( Kowalsk i et al.|2008 ) as imple- 



mented in CosmoMCQ ( |Lewis fe Bridle| |2002 



Lewis & Bri- 



|dle| with systematic errors flag turned on, and BAO to the 
two data points r s (zd ec )/ Dy{z — 0.2) and r s (zdec)/ Dy(z = 
0.35) from Percival et al. (2007). For reference we have done 
also 6 MCMC runs (excluding WMAP&ACBAR from the 
above list) for the spatially flat non-interacting model. Each 
of these 13 runs has 3-4 chains with mean input multiplic- 
ity in the range 3-10, and the number of accepted models 
in each chain is ^25000. The measure of mixing, the worst 
eigenvalue R — 1 (which is better the closer it is to zero), is 
for all cases less than 0.03. More technical details are given 
in Appendix [A] 

Table [A2] on page [14] in Appendix [B] collects the best- 
fitting models of each run (excluding WMAP& ACBAR). 
In addition, we show for each data-set two relatively 
good- fit models with a strong negative/positive interaction 
(|r/iJo| > 0.1) which are not 'far' from the best-fitting non- 
interacting model - in a sense that A% 2 = X 2 (strong inter- 
action model) - x 2 (best-fitting model) < 4. In the last four 
rows of the table (Models 25-28) we show interacting models 
which have most of their parameters equal to the WMAP 
best-fitting non-interacting model (Model 1). 



http:/ /cosmologist.info/cosmomc 
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Figure 2. Models with most of their parameters equal to the parameters of the best-fitting to WMAP non-interacting model. 



Fig.[2]shows the angular power spectra, distance modu- 
lus, and r s (zdec)/ 'Dv(z) for Models 25-28 and Model 1 from 
Table [A2] The WMAP best-fitting non-interacting model 
(Model 1; thin green/grey line) fits well all the data. The 
only exceptions are that it fails to fit the low quadrupole 
I — 2 in the WMAP TT spectrum, and undershoots both 
the z — 0.20 and z — 0.35 BAO data points, whose error bars 
are smaller than the asterisk symbol in the plot. Therefore 
it would be surprising if the interacting model could fit the 
data overall much better. Now we take Model 1 and turn the 
interaction on to |r/fZo| = 0.3 (Models 25 and 26). The neg- 
ative (positive) interaction leads to an extremely bad fit to 
the CMB TT data, as the the model now vastly overshoots 
(undershoots) the first and second acoustic peaks. This is 
also reflected in the TE spectrum where the peak at I ~ 300 
is overshot (undershot). The SN data cannot distinguish the 
interacting model from the non-interacting model, though 
there is a tiny difference: a negative (positive) interaction 
leads to a faster (slower) increase of the distance modulus 
as a function of redshift than seen in the non-interacting 
model, and hence to a slightly better (worse) fit to the SN 
data. Finally, a negative interaction helps fitting the BAO 
data as it systematically shifts r s (zdec)/ 'Dv(z) upward. A 
positive interaction makes the fit to BAO data worse. All 



these remarks are reflected in the \ 2 values presented in 
Table lA2l 

It is quite straightforward to understand why a nega- 
tive interaction helps in fitting the SN and BAO data. The 
negative interaction means that the energy transfer is from 
dark energy to dark matter. As we have kept today's val- 
ues of co c (and Qdeo) and Ho fixed, this means that in the 
negatively interacting model there has been more dark en- 
ergy in the past than in the non-interacting model. This 
causes more acceleration of the expansion of the universe, 
and hence larger distance moduli at high redshifts. In the 
BAO the distance measure is in the denominator, and at 
first sight one would expect smaller r s (zd ec )/ 'Dv(z) than in 
the non-interacting case. However, as Dy is proportional to 
L>y s the effect of an increased angular diameter distance, 
Da(z), is mild. Instead the main effect now comes from the 
different sound horizon. For a fixed today's u c the negatively 
interacting model has a smaller cold dark matter density at 
last scattering uj c (z*) than the non-interacting model. The 
sound horizon r s (zd ec ) is proportional to where /c eq is the 
wave number corresponding to the scale A eq that re-enters 
the horizon at matter-radiation equality. The smaller uj c (z) 
at early times means that the matter density exceeds the ra- 
diation energy density later than in the non-interacting case. 
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As small scales re-enter the Hubble horizon before the large 
scales, A eq (or k~^) and hence the sound horizon r s (zd ec ) 



will be larger in the negatively interacting case. Table |A2| 
confirms these conclusions (compare Dy and r s in Models 1 
and 25). The positive interaction model behaves vice versa. 

Now we can explain also the bad fit to the CMB. The 
angular power spectrum is very sensitive to the dark mat- 
ter density at last scattering oj c (z*) (or at the redshift of 
matter-radiation equality z eq ). As the negative interaction 
model has a small uj c (z*) and hence a small z eq , the matter- 
radiation equality occurs very close to last scattering z*. 
Therefore last scattering does not happen in the matter 
dominated era but around the transition from radiation 
domination to matter domination. This causes a large early 
integrated Sachs- Wolfe effect (elSW), which amplifies the 
first (and second) acoustic peaks. The angular power spec- 
trum looks like that of a non-interacting model with a very 
small today's uj c . The positive interaction model looks like a 
non- interacting model with a very high today's cj c , i.e., the 
first and second acoustic peaks are vastly undershot (due to 
the early matter-radiation equality and hence a small elSW 
amplification). Therefore, in order to obtain a good fit to 
the CMB in the interacting model, we need to adjust to- 
day's uj c in such a way that cu c (z*), or rather a eq , matches 
with the best-fitting non-interacting model. We have done 
this in Models 27 and 28; see Fig. [2] and Table [A2] In the 
best-fitting non-interacting model (Model 1) today's dark 
matter density is cj c = 0.107, whereas to obtain a good fit 
with a negative (positive) interaction of T/H = —0.3 (+0.3) 
we require a larger oo c = 0.137 (smaller oo c = 0.080). As 
the larger (smaller) dark matter density leads to a smaller 
(larger) Qdeo, the interacting model will now have a smaller 
(larger) angular diameter distance to last scattering Da(z*). 
This would lead to all the acoustic peak structure shifting 
slightly to the left (right) from the WMAP data. As D A 
is proportional to Hq 1 , this mismatch can be corrected by 
decreasing (increasing) the value of Ho so that we obtain 
roughly the same Da(z*) and hence the same acoustic peak 
positions as in the non-interacting case. With a negative 
(positive) interaction of Y/Ho = —0.3 (+0.3), we require 
H = 56 (66), while the best-fitting non- interacting m odel 
has Ho = 61kms _1 Mpc _1 . (Note: D v and r s in Table |A2 
are reported in units /i _1 Mpc. The interacting Models 27 
and 28 have r s — 154 Mpc, which is exactly the same as for 
the non-interacting Model 1.) 

Fig.[2]and Table [A2] show that after the above-described 
adjustments the strongly interacting \T/Hq\ — 0.3 Models 27 
and 28 provide an excellent fit to the CMB acoustic peaks. 
In the high-/ (I > 32) region of the TT and all of the TE 
spectrum, the interacting models are totally indistinguish- 
able from the non-interacting best-fitting Model 1. However, 
as the negatively interacting Model 27 has small Qdeo = 0.49 
and hence also a smaller Qde in the recent past than the non- 
interacting model, the late integrated Sachs-Wolfe effect is 
suppressed, so there is less TT power at low multipoles. As 
the non-interacting model slightly over shoots the Cj T spec- 
trum at low multipoles, the negatively interacting model 
leads to a better fit her e; se e Fig. [2] upper left panel and the 



y 2TT 
Z<32 



in Table 



A2 



for Models 1 and 27. (Note: the 



column Xi< 

Xi<32 numbers are from Gibbs sampling of the actual CMB 
map, and the more negative the number is the better the 
fit is.) The positive interaction Model 28 has a larger Qde in 



the recent times than the non-interacting model, and hence 
gives rise to a large late ISW effect, and a poor fit to low 
multipoles in the TT spectrum. 

The SN data cover the relatively recent past only and 
hence the small Qdeo = 0.49 in the negative interaction 
Model 27 leads to a poorer fit. The SN data indeed slightly 
favour more acceleration in the recent past and the posi- 
tive interaction Model 28 (^d e o = 0.76) has this property, 
leading to a better fit to SN than the non-interacting Model 
1, which has Qdeo — 0.6 5; se e Fig. [2] lower left panel and 



the column Xsn m Table A2 For the same reason the neg- 



ative interaction model fits the BAO data worse than the 
positive interaction model. Also the BAO data cover only 
the relatively recent past and favour more acceleration be- 
tween redshifts z = 0.20 and z = 0.35 than the best-fitting 
non-interacting model provides. 

Indeed already from Table |A2] we can see these general 
trends. All the best-fitting interacting models (2, 14, 18, 
and 22) have negative F and are slightly better fits to the 
CMB due to sightly better fit to low-/ CMB spectra. In gen- 
eral the best-fitting interacting models have indistinguish- 
able xIn a nd Xbao from the best-fitting non-interacting 
models. The largest- |T| models within A\ 2 < 4 from the 
best-fitting model (see Models 3, 4, 15, 16, 19, 20, 23, 24 
in Table A2) always obey the following: a good- fit large 
negative interaction model is by about A% 2 = 1-3 better fit 
to WMAP than a good-fit large positive interaction model, 
again due to the low-/ behaviour. However, a negative in- 
teraction model is always a worse fit to both SN and BAO 
data (due to too small Qde in the recent past). 

The ISW effect in the temperature angular power spec- 
trum is given by 



ISW 



xi ( (fc(r -r))e- T ° d(r ' T0) ]drJ ^ , (6) 



where ji is the spherical Bessel function and T Q d(T, To) = 
JJ° opacity (r)df is the optical depth from r to today (to). 
Note that at late times (0 < z < 10) we have $ = — ^, 
since the anisotropic stress vanishes. As it is the combina- 
tion ^' — which defines the ISW effect, we show on the 
left panel of Fig.|3]the evolution of ^ — <I> for the same mod- 
els as in Fig.|2j i.e., Models 1 and 25-28 from Table |A2] At 
high redshifts z > 10 4 , deep in the radiation era, the evolu- 
tion in the interacting models is indistinguishable from the 
non-interacting model. If we keep all the other parameters 
fixed to the best-fitting non-interacting model, then a neg- 
ative (positive) interaction leads to faster (slower) decay of 
the potential ^ — $ around the time of last scattering; com- 
pare the dashed blue (solid black) curves to the green curve 
on the left panel of Fig. [3] This matches to what we already 
explained about the early ISW effect: it is more (less) pro- 
nounced in the case of negative (positive) interaction, since 
the matter-radiation equality appears later (earlier) than in 
the non-interacting case. At late times, ^ z < 10, the in- 
teraction starts to modify the ISW source directly. The neg- 
ative (positive) interaction leads to more gradual (steeper) 
decay of ^ — and to a smaller (an enhanced) late ISW 
effect. 

If we adjust oo c and H (and hence Qdeo) so that we ob- 



6 Jussi Valiviita, Roy Maartens and Elisabetta Majerotto 



1.3 



1.2 



1.1 



0.9 



0.8 



0.7 




Model 25 (As Model 1 , but r/H Q = -0.3) 

Model 26 (As Model 1 , but r/H Q = +0.3) 

. _ . _ . Model 27 (as 25, but co c and H Q adjusted) 

Model 28 (as 26, but co c and H Q adjusted) 
- Model 1 (best-fitting WMAP non-inter.) 



ZE 



ZE 



10 



10 
z 



10 



1.3 



1.2 



1.1 



i 



0.9 



0.8 



0.7 




Model 27 
Model 28 

As Model 27, but r=0 in perturbation eqns 
As Model 28, but r=0 in perturbation eqns 
Model 1 (best-fitting WMAP non-inter.) 



ZSZ 



10 



10 
z 



10 



Figure 3. Redshift evolution of the ISW source. Dashed vertical lines indicate last scattering (z 
interaction starts to directly modify the evolution of the ISW source (z « 10). 
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tain the same z eq as in the non-interacting case, and hence 
a perfect fit to the acoustic peaks in the data, the above- 
described effects on the late ISW effect become even more 
pronounced: compare the dot-dashed red (solid cyan) curves 
to the green curve on the left panel of Fig. [3] This is because 
with negative (positive) F we need a smaller (larger) Qdeo, 
and so the background effect on the gravitational poten- 
tial is to reduce (increase) its decay rate. Now an interest- 
ing question arises: how much of the late ISW effect in the 
good-fit models comes from the different background, i.e. 
different Qdeo compared to the non-interacting best-fitting 
model, and how much comes from the modified perturbation 
evolution equations. The right panel of Fig. [3] addresses this 
question. Blue dashed (black solid) curves show how the neg- 
ative interaction Model 27 (positive interaction Model 28) 
would behave if we ignored the interaction in the pertur- 
bation equations. Interestingly, we would drastically over- 
estimate the effect of interaction on the late ISW effect. In 
the perturbation equations, see equations (11-14) or (29-41) 
in the companion paper (Majerotto, Valiviita & Maartensl 



2009), there seems to be a term which partially cancels the 



effect of different background evolution. This is easiest to 
see in the longitudinal (conformal Newtonian) gauge where 
B — = E. As the dark energy perturbations remain sub- 
dominant, the interaction terms in their evolution equations 
cannot be responsible for the difference. Moreover, in the 
CDM velocity equation explicit interaction term is com- 
pletely missing. So the only perturbation equation where 
the interaction appears to have an effect is the CDM den- 



sity contrast equation, which in the longitudinal gauge reads 

6' c = 3?// - aYct) . (7) 

As long as the scale factor is small (a < 1) the evolution 
is like the non-interacting evolution, apart from the differ- 
ent evolution of the background, but in the later matter 
dominated era or dark energy dominated era, when a starts 
to approach 1, the interaction starts to modify directly the 
way the CDM reflects gravitational potential wells. In this 
era there is no anisotropic stress, so <j) — ip, and let us as- 
sume that ip > 0. Then due to the background effect of dark 
energy, ip' < 0. Now the effect of the — aT<j> term is that if 
r < (> 0), then the CDM density contrast decays slower 
(faster) than in the non- interacting case. This in turn feeds 
back into the evolution of via an Einstein equation (the 
general relativistic Poisson equation, equation (39) in the 
companion paper (Maj erotto, Valiviita &; Maartens||2009 )). 
The effect is stronger the closer to 1 the scale factor is. 

In the Newtonian gauge we can always neglect the 
dark energy perturbations, e.g. the total density contrast is 
5 = Spc/ptot at late times. However, in synchronous gauge 
we cannot do this for the velocity perturbation since C = 0. 
Indeed, at late times to t ~ Ode- But the second Einstein 
equation is k 2 r)' s = AitGa 2 (p + p)0t o t - see e.g. equation 



(21b) in Ma & Bertschinger (1995) - and this is how the 
synchronous gauge metric perturbation rj s is calculated in 
CAMB ( |Lewis et al-pOOO] ). So in synchronous gauge the in- 
teraction in the perturbation equations enters the ISW effect 
via the interaction in the 0' de equation, (12) in the compan- 
ion paper ( |Majerotto, Valiviita &; Maarte ns 2009). 

Note that it is difficult to go to very large positive in- 
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Figure 4. 2d marginalized likelihoods for the interacting model with only the WMAP, only the SN, and only the BAO data. The darker 
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differs from the tabulated one which is at T/Hq = 2.92; see Table |A2| 
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Figure 5. 2d marginalized likelihoods for the interacting model with the WMAP, WMAP&SN, WMAP&BAO and WMAP&SN&BAO 
data. 



8 Jussi Valiviita, Roy Maartens and Elisabetta Majerotto 







13.5 



14 14.5 15 

Age/GYr 



WMAP only; ■ best-fit 

WMAP&SN; • best-fit 

WMAP & BAO; ▼ best-fit 

WMAP & SN & BAO; ♦ best-fit 



Figure 6. Id marginalized likelihoods for the interacting model with the WMAP, WMAP&SN, WMAP&BAO and WMAP&SN&BAO 
data. 



teractions, since the larger the T, the larger w e one needs 
in order to avoid negative Qde in the past, i.e., to avoid 
the zero crossing of pde which causes the perturbation equa- 
tions becoming singular at that moment, as discussed in 
Appendix |A3| So we expect the large interaction models to 
fit the CMB data as well as the non-interacting model, with 
negative Y sightly favoured since this improves the fit in the 
low multipole region. Moreover, in marginalized likelihoods 
the negative interaction models will be favoured, since the 
pde zero-crossing problem means that there is much less vol- 
ume in the allowed parameter space in the positive Y region 
than in the negative Y region; see Appendix |A3| There will 
be a strong degeneracy between cj c and T, Ho and T, and 
Qdeo and Y. When adding the other data (SN or BAO or 
SN&BAO) we expect most of the negative interaction mod- 
els that fit the CMB alone well, to be excluded due to their 
very small Qde today and in the recent past. 



3 LIKELIHOODS 

As predicted, due to the degeneracy between u c and 
Y/Ho, the CMB data alone do not provide tight con- 
straints on the interaction. The degeneracy is almost lin- 
ear, (jO c ~ 0.107 - O.ir/ifo, according to Fig. [I] which 
shows 68% and 95% confidence level (CL) regions with 
WMAP, WMAP&ACBAR, and WMAP&SN&BAO data. 
As explained in the previous section, once cj c (and Ho) are 
adjusted, the interacting model produces completely indis- 
tinguishable CMB angular power spectra at I > 32. There- 
fore, combining WMAP (which reaches up to / ~ 1000) with 
ACBAR (which reaches / ~ 2000) does not help at all. In- 
deed, even with the forthcoming Planck data one will not be 
able to improve the constraints presented in Fig. [I] unless 
supplemented with some other non-CMB data. This is be- 
cause the only signature from the interaction appears in the 
ISW region, and there the accuracy of CMB data is already 
now cosmic variance limited. 

Furthermore, we see from Fig. ^ that according to the 
marginalized likelihoods with the CMB data, negative inter- 
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Figure 7. Comparison of Id marginalized likelihoods for the interacting and non-interacting models with the WMAP and WMAP&SN 
data. The solid lines are for the non-interacting (T = 0) reference model. 



actions are strongly "favoured" over positive interactions. 
As explained in Appendix |A3| this is partially due to a vol- 
ume effect caused by shrinking of the allowed w e direction 
of parameter space for large positive interaction. In addi- 
tion, positive interactions worsen the fit to the ISW region, 
while negative interactions improve it. According to Fig. ^ 
the 95% CL region of T/H extends from -0.9 to +0.1. How- 
ever, it should be noted that even stronger than -0.9 negative 
interactions would become allowed, if we lowered the lower 
bounds of two of our top-hat priors: 40 < H < 100, and 
< Qdeo < 1. This becomes evident later in the second and 
third panels of Fig. [5] 

Adding SN and BAO data to the analysis leads to a 
more symmetric 95% CL region -0.23 < T/H Q < +0.15, 
as seen in Fig. ^ In this case, the worse fit to ISW and 
the volume effect from w e in the case of positive interaction 
become cancelled by better fits to the SN and BAO data. In 
order to gain more insight into this, in some cases dangerous, 
competition between the CMB and other data, we compare 



in Fig. [4] selected 2d marginalized posterior likelihoods when 
using only WMAP or only SN or only BAO data. 

The SN or BAO data alone do not significantly con- 
strain any parameters of our model other than those shown 
in Fig. [4] i.e., they push oo c down to oo c < 0.15, Qdeo up to 
0.65 < tt d eo < 0.90, and w down to w < -0.75 at 68% 
CL. The non-interacting model (r = 0) is consistent with 
all three data sets (WMAP, SN, BAO). Most interestingly, 
the non-interacting model sits in the intersection of all three 
data sets so that there is no tension between them. As al- 
ready noticed, there is a tension between CMB and SN or 
CMB and BAO in the negatively interacting models. This 
tension is most pronounced in the top right panel (r, Qdeo) 
of Fig. [4] However, the situation is not too bad since there 
is plenty of parameter space volume in the intersection of 
95% CL regions of WMAP and SN or WMAP and BAO. 
Although the SN or BAO data do not put any direct con- 
straints on the interaction (even |T/iJo| 00 3 fits them well) 
the net effect of combining CMB with SN or BAO data is 
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to force Qdeo > 0.65 and hence to cut away the large neg- 
ative interaction models, leaving an almost symmetric re- 
gion around r = 0. This effect is seen in the top left panel 
of Fig. [5] where we show 2d marginalized likelihoods from 
our MCMC runs for the interacting model with the WMAP 
data alone, the WMAP&SN, the WMAP&BAO, and the 
WMAP&SN&BAO data. The SN and BAO data are rather 
consistent with each other when constraining the interacting 
model. Therefore combining WMAP with either or both of 
them leads to very similar constraints, as seen in Fig. [5] 

In Fig. [6] we show the Id marginalized likelihoods for 
all of the primary MCMC parameters of our model and for 
two derived parameters: Qdeo and the age of the universe. 
We note that WMAP combined with the BAO data prefers 
slightly larger today's Hubble parameter Ho and dark energy 
density Qdeo than with the SN data. The CMB data alone 
favour negative interactions and thus require small H and 
Qdeo- This corresponds to a very old universe as seen in the 
last panel of Fig. [6] 

As mentioned above, from the Id plot for V/Hq in 
Fig. [6] the negative interaction seems more probable than 
the positive interaction. Indeed with the WMAP data 
(WMAP&SN&BAO data) 96.8% (77.6%) of models in 
our Markov Chains have a negative T, which means en- 
ergy transfer from dark energy to dark matter. However, 
we should be cautious in claiming that the WMAP data 
(WMAP&SN&BAO data) favour energy transfer from dark 
energy to CDM at 96.8% CL (77.6% CL). It should again be 
stressed that this is partially the volume effect from the w e 
direction of the parameter space. Indeed, recently [Pereira fe| 
Jesus| ( |2009| claimed tha t with 93% probability the data - 
which in |Pereira &; Jesus ( 2009 ) were the background-based 
data only - favour decay of dark matter to dark energy. 
Now we have about the same probability in favour of energy 
transfer from the dark energy to dark matter. The message 
here is that, in addition to the volume effect, these claims 



are highly model dependent: Pereira & Jesus (2009) studied 



an interaction proportional to pde whereas our interaction is 
proportional to p c . 

We summarize here our most stringent results for the 
interacting model by giving minimal 95% intervals ( Hamann 



eTaLl[2007t from our MCMC run with WMAP&SN&BAO 
data: u b e (0.0212, 0.0241), u c G (0.859, 0.125), H € 
(63, 70), r G (0.057, 0.133), F/H G (-0.23, +0.15), w e G 
(-0.80, -0.19), w G (-1.00, -0.63), n s G (0.937, 1.002), 
ln(10 10 A|) G (2.95,3.14), Q de o G (0.648,0.767), AgeG 
(13.6, 14.3) Gyr. 

In Fig. [7] we compare the Id marginalized likelihoods 
of the interacting model to the non-interacting reference 
model (r = 0). The key differences are: the interacting 
model leads to broader distributions of lj c , Ho and Qdeo- 
This is due to the degeneracy between u c (or Ho or Qdeo) 
and r. Moreover, smaller H is favoured by the interacting 
model. However, note again that this is partially a parame- 
ter space volume effect, since in our MCMC chains there are 
many more negative interaction models than positive inter- 
action models, and the good- fit negative interaction models 
have small Ho- Finally we note that Ho, both in the in- 
teracting and non-interacting cases, is smaller than in the 
ACDM model, where a typical result of a likelihood scan 
peaks around Ho = 72 kms _1 Mpc _1 . This is because we let 
wo (and w a ) vary, and there is a strong degeneracy (even 




Figure 8. Comparison of a 2d marginalized likelihood for the 
interacting and non-interacting models with the WMAP and 
WMAP&SN data. Note the degeneracy between Ho and wo- The 
solid lines are for the non-interacting (r = 0) reference model. 



in the non-interacting case) between Ho and wq (or Wde)- 
From Fig. [8] we see that a wo = —1 model prefers largest Ho 
(and if we allowed wo < — 1, even larger values of Ho would 
be favoured). 



4 ADDITIONAL CONSTRAINTS ON THE 
MODEL 

In the two previous sections we have seen that any CMB 
data alone fail to constrain the interacting model, in partic- 
ular with negative interaction, but the constraints from the 
SN or BAO data for the background dark energy density in 
the recent past cut out the strong interaction models with 
r < 0, as these models have very small Qde- In this section 
we briefly discuss some other data that could be used. 

Most interestingly, as the interaction seems to affect 
the ISW effect only, accurate CMB - large-scale structure 
(LSS) cross-correlation data on the ISW effect (combined 
with other data sets such as CMB, SN and BAO) may turn 
out to set the most stringent constraints on the interaction. 
Firstly, the CMB-LSS cross-correlation data probe the evo- 
lution of Qde over time in the recent past. Secondly, as the 
interaction affects not only the background evolution, but 
directly the evolution of perturbations (in particular, the 
way that the evolution of the cold dark matter perturba- 
tion and gravitational potential are linked to each other) 
at redshifts probed by the ISW data, the ISW data may 
see the effects from the interaction most directly. We are 
working on this issue. While we were finalizing this paper, a 
work on other type of dark sector interaction appeared also 
suggesting that the interaction could be detected via its ef- 
fects on the ISW signal (He et al. 2009). Earlier Lee et al 
(2006} discussed the modified ISW effect in an interacting 
quintessence model. However, they fixed all the other pa- 
rameters to (or near to) the best-fitting to WMAP Lambda 
CDM model. Therefore, although mentioning the degener- 
acy between Qdeo (or uj c ) and the interaction rate elsewhere 
in their paper, they missed the degeneracy when finding con- 
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straints on the interaction rate with WMAP data, ending 
up with artificially tight constraints on the interaction with 
CMB data alone. This can be seen, for example, in the third 
panel of our Fig. [5] Had we fixed Qdeo to 0.76, as |Lee et al.| 
(2006) did, our constraint on the interaction rate would be 
misleadingly tight: -0.03 < T/H < +0.005. It should be 
noted that although adding SN or BAO or SN&BAO into 
the analysis forces Qdeo (and oo c and Ho) closer to the best- 
fitting (to WMAP) Lambda CDM model and hence leads 
to a similar kind of effect as fixing the parameters, letting 
the other parameters freely vary and consequently taking 
properly into account the degeneracy gives much looser con- 
straint: -0.23 < T/Ho < +0.15 with WMAP&SN&BAO. 
So fixing the other parameters would misleadingly give more 
than an order of magnitude smaller 95% interval in our case. 

Other interesting data come form the galaxy-galaxy 
power spectrum which is a probe of the underlying theoret- 
ical matter power spectrum. However, we have decided not 
to use these data here, since the exact relation between the 
observed galaxy-galaxy power spectrum and matter power 
spectrum is not known (due to redshift space distortions, 
non-linearities, etc) and moreover we calculate the pertur- 
bation evolution in the linear regime only. Probably the in- 
teraction would affect non- linear structure formation (Baldi 
et al. 2008; Baldi 2009). So at the moment it does not make 
sense to fit the shape of the (linear) matter power to the 
data. For similar reasons we do not consider weak-lensing 
data. 

As the strongly interacting Models 27 and 28 in Ta- 
ble |A2| have very different recent-time evolution of matter 
and dark energy densities, the overall amplitude of the mat- 
ter power is affected. This results in different ag (root mean 
square mass fluctuation on 8/i _1 Mpc scale) than in the non- 
interacting model. We find that today as(z = 0) = 0.598, 
0.628, 0.654 for Model 27 (F/H = -0.3), Model 1 (r = 0), 
and Model 28 (T/Ho = +0.3), respectively. Therefore, we 
predict that as measurements could rule out large negative 
interactions as they may lead to too small erg. On the left 



panel of Fig. [5] we show the redshift evolution of as- We 
have normalized the curves for the interacting models to 
fit today's ag of the non-interacting model. For example, 
for the T/Ho — —0.3 model as(z = 0) is smaller than in 
the non-interacting case, and towards the past it decreases 
even faster than in the non-interacting case. On the right 
panel of Fig. [5] we compare matter power spectra at red- 
shifts z — and z — 2 in the strongly interacting and non- 
interacting models. Small-scale (large k) shape is unaffected 
by the interaction, whereas large-scale (small k) shape (and 
amplitude) are affected significantly. Marginalizing analyt- 
ically over the galaxy bias, the SDSS DR4 Luminous Red 
Galaxies sample ( |Tegmark et al. 2004| |2006 ) gives the best 
X 2 for the positive interaction and worst x f° r the negative 
interaction, while the x 2 of the non-interacting model falls 
between these two. As the effect of adding the LSS data to 
the analysis seems to be very similar to the effect of adding 
BAO or SN, we expect that adding the LSS data would give 
more weight to the positive interactions and hence lead to 
even more symmetric and tighter probability distribution 
about r = than with the CMB&SN&BAO data. 



5 CONCLUSION 



In the companion paper (Majerotto, Valiviita & Maartens 



2009) we have presented, for the first time, a systematic 



derivation of initial conditions for perturbations in inter- 
acting dark matter - dark energy fluid models deep in the 
radiation dominated era. These initial conditions are essen- 
tial for studying the further evolution of perturbations up 
to today's observables. We have focused on the interaction 
Q% = — rp c (l + 5 c )u%, where T is a constant which has the 
same dimension as the Hubble parameter H] see Eqs. |l]) and 
(g). In our previous work (Vali viita et al.||2008| we showed 
that if the equation of state parameter for dark energy is 
— 1 < Wde < —4/5 in the radiation or matter dominated 
era, the model suffers from a serious non-adiabatic insta- 
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bility. However, in this paper (and in the companion paper 
(Majerotto, Valiviita & Maartens 2009 )) we have shown that 



the instability can easily be avoided, if we allow for suitably 
time- varying dark energy equation of state. Our worked out 
example is for the parametrization Wde = woa + w e (l — a). 
With this parametrization, viable cosmologies (in the inter- 
acting model) result whenever wo is close to — 1 and w e < — 1 
or — 4/5 ^ w e < 1/3, as long as wq + 1 and w e + 1 have the 
same sign. 

We have implemented into a publicly available Boltz- 
mann integrator, CAMB, the background equations, and 
first order (linearized) perturbation equations, as well as 
the adiabatic initial conditions for the interacting model 
with time- varying equation of state parameter. We have 
performed full Monte Carlo Markov Chain likelihood scans 
for this model as well as for the non-interacting (r=0) 
model for a reference, with various combinations of publicly 
available data sets (WMAP, WMAP&ACBAR, SN, BAO, 
WMAP&SN, WMAP&BAO, WMAP&SN&BAO). To avoid 
ending up with too complicated shape of parameter space 
we have focused on non-phantom models, —1 < Wde < 1/3. 
(In addition, we consider phantom models to be unphysical.) 

The main result is that there is a degeneracy between 
the interaction rate and today's dark energy (or dark mat- 
ter) density in light of CMB data. Therefore, CMB data 
alone cannot rule out large interaction rates, not even 
Planck, since the high-multipole part of the CMB angular 
power spectra are totally indistinguishable from the non- 
interacting case. The only signal remaining from a large in- 
teraction rate would be a modified integrated Sachs- Wolfe 
effect, which makes using the ISW data an appealing line 
of future work. In this paper we have broken the degener- 
acy by Supernovae data and by baryon acoustic oscillation 
data, finding that the CMB&SN&BAO data constrain the 
interaction rate to about 20% of the expansion rate of the 
Universe. 
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APPENDIX A: TECHNICAL DETAILS 
Al The code and modified sound horizon 



We have modified publicly available CosmoMC (Lewis & 
Bridle|2002[|Lewis fe Bridle] and CAMB ( |Lewis et al.|2000 f 
for this study. Into CAMB we have implemented the in- 
teracting evolution equations for the background and 
([2]), as well as the interacting perturbation evolution for 
dark energy, equations (11-14) from the companion paper 
(Maj erotto, Valiviita &; Maartens|2009 ), in the synchronous 
gauge: B — 4> — 0, ip — rj s , and E — — k~ 2 (6r] s + h s )/2, with 
7] s and h s representing the synchronous gauge metric per- 
turbations. In the synchronous gauge, the perturbed CDM 
equations of motion appear to look the same in the inter- 
acting and non-interacting cases, since = 0. However, it 
should be noted that as the background evolves differently, 



it affects the evolution of h s , and hence S c . We have imple- 
mented the initial conditions for perturbations deep in the 
radiation dominated era, specified in the companion paper 



(iMajerotto, Valiviita & Maartens||2009 ). 



In order to use the BAO data in the interacting model 
an additional modification to CAMB is necessary. In the 
standard version of CAMB/CosmoMC, the sound horizon is 
calculated by numerically integrating up to last scattering 



PT* pa* 

s (a*) = / c s (r)dr — / c s (a)/ada, 
Jo Jo 



(Al) 



where c s is the sound speed in the photon-baryon fluid and 
a* is very accurately calculated from a fitting formula (Hu 



&; Sugiyama| |l996) which is valid if the matter scalings are 
the standard ones, pb oc a -3 and p c oc a -3 . Firstly we need 
a different formula, as for BAO we want to integrate up to 
decoupling adec, for which one has another fitting formula 
"T996J) 



(Hu & Sugiyama 



Also this formula is valid only if 
However, we know that this is not 
adec numerically, and 



pb oc a " and p c oc a 
true for our case. Therefore we find 
then calculate numerically 

r.s(aaee) = / c s (a)/a da 



f 

Jo 



(A2) 



The defining equation of decoupling is (Hu & Sugiyama 



1996) 



/' 

J Tn 



R(r) x opacity(r) dr = 1 . 



(A3) 



where R = ^pb/pj- We numerically follow this integral from 
today toward past times f until the value 1 is reached. Then 
we record the value of the scale factor at this moment, name 
it adec and convert to a redshift Zdec = 1 / ^dec — 1 • Finally, to 
match the definitions in Eisenste in et al.| ( |2007| and |Percival| 
et al.| ( [2007| (from where we take the BAO data) we multiply 
the result by their "phenomenological" factor 0.96; Zdec = 
0.965dec- We have verified that in the non-interacting case 
these definitions and our numerical routines lead to the same 
and r s (adec) as those given as a test case in Percival 



Zd 



eTaT1 ( [2007| 



A2 The parameters and their prior ranges 

In the interacting model we have 9 primary MCMC param- 
eters which we vary over wide ranges with uniform (flat) 
prior over their range: the physical baryon density today 
ujb — h 2 Qbo, the physical cold dark matter density today 
co c = /i 2 Q c0 , the Hubble parameter today ifo, optical depth 
to reionization r re , the interaction F in units of today's Hub- 
ble parameter, the early dark energy equation of state pa- 
rameter w e , the dark energy equation of state parameter 
today wq, scalar spectral index ns of the primordial pertur- 
bations, and the amplitude of primordial perturbations As 
[we use log e (10 10 A|) as in the standard CosmoMC]. In the 
non-interacting reference model we have 8 free parameters 
as we keep T fixed to zero. 

We exclude phantom models (wde < — 1) as unphysical, 
and we focus on the following ranges 

-l<w o <0, -l<w e <0.3. (A4) 

It should be noted that the CMB data actually exclude all 
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Figure Al. The K; e volume effect. The circles show samples from our MCMC run with the WMAP&BAO data. The colour scale 
indicates the value of wq for each sample. 



blow-up models, i.e., the interacting models with — 1 < w e < 
—0.8. Keeping this in mind, our motivation to drop phan- 
tom models from the analysis becomes reinforced: it would 
make the posterior of w e consist of two separate intervals, 
and hence the posterior likelihood of w e would have two 
peaks, one with w e < — 1 and another one somewhere in the 
range —4/5 < w e < 1/3. As the MCMC technique is par- 
ticularly ill-suited for multiple peaked posteriors, allowing 
for phantom models would ruin the analysis. Moreover, we 
should demand w e + 1 and wo + 1 to have the same sign, 
since otherwise there would be a time r oc in the past where 
Wde(aoc) = — 1. Such a 'minus one' crossing causes severe 
problems since the dark energy perturbation equations con- 
tain terms proportional to 1/(1 + Wde), see the companion 
paper ( |Majerotto, Valiviita &; Maart ens 2009). The require- 
ment signfwe + 1) = sign(i(;o + 1) would further complicate 
the shape of the parameter space, if we allowed for w e < — 1 
or wo < —1. 

As discussed in | Valiviita et ah (2008), 
would conflict with the data. We find that 
range is T/Ho G (—4, +4). The ranges of other remaining 
primary MCMC parameters are: cob G (0.005, 0.1), oo c £ 
(0.01, 0.99), Ho e (40, 100), r e (0.01, 0.4), n s G (0.5, 1.5), 
ln(10 10 A|) G (2.7, 4.0). Finally, we have an amplitude of 
the SZ template in CosmoMC, Asz G (0, 2), which is used 



\T/H \ > 1 
a reasonable 



when adding the SZ templates to the high-Z part of the TT 
angular power spectrum. We also restrict the analysis to 
positive dark energy densities today, i.e., apply a top- hat 
prior Qdeo > 0. As we study spatially flat models, the mat- 
ter density today is Q m o = 1 — ^d e o- 

In Table |A1| we list some of the above-mentioned pa- 
rameters as well as other symbols used in this paper. 



A3 Positive V and zero-crossing of pa e 

There is one additional complication in studying the inter- 
acting model and interpreting the marginalized posterior 
likelihoods. As pointed out in Valiviita et al. ( 2008| , in the 
case of constant Wde, any positive interaction T/Ho > 
would lead to a zero crossing of pde, or in other words Qde- 
This means that starting the background calculation from 
today's positive value, say Qdeo ~ 0.7, and integrating back- 
ward in time, at some moment Td e ,zc in the the past Qde(r) 
crosses zero and becomes negative for r < Td e , zc - While we 
lack deep understanding of the nature of dark energy we 
might even accept this possibility. However, the perturba- 
tion equations (11) and (12) in the companion paper (Ma- 



jerotto, Valiviita & Maartens 2009), have the dark energy 



density in the denominator, and therefore become singular 
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symbol explanation denning equation / reference 



subscript c cold dark matter 

subscript de dark energy 

subscript b baryons 

subscript variable evaluated today 

r conformal time dr = a~ 1 dt, where a is the scale factor of the Universe 

' conformal time derivative E.g. a' = ^ 

H Hubble parameter, i.e., the expansion rate of the Universe H = ^ /a 

TL conformal Hubble parameter Tt = a' /a = aH 

uob physical density parameter of baryons uo^ = flboh 2 

uj c physical density parameter of cold dark matter uj c = Q c oh 2 

ptot total density sum of energy densities of all constituents of the Univ. 

ns scalar spectral index 

As amplitude of the primordial curvature perturbation power V-jz(k) = As(/c//co) ns_1 , where ko = 0.05 Mpc -1 

r re optical depth to reionization 

Dy dilation scale equation (2) of Eisenstein et al. (2005) 

Da angular diameter distance 

Zdec redshift of decoupling Eq. (A3) and explanation after that: z^ ec = 0.965dec 
redshift of last scattering 

r s sound horizon a decoupling Eq. ( A2 ) 

4>, tpj B, E scalar metric perturbation variables equation (6) of iMajerott o et all |2009| 

<E>, ^ scalar gauge invariant metric perturbations equations (6) and (25) of Majerot to et al.N2009j 

?7 S , h s synchronous gauge metric perturbations these ar e called 77 and h inTO a & Bertschinger (1995) 

6 C cold dark matter velocity perturbation see e.g. ( | Valiviita et al.| 2008 Ma & Bertschinge r |l995| ) 

Ode dark energy velocity perturbation 



Table Al. List of selected symbols used in this paper, with their physical meaning and defining equation or a reference. 



at the moment Td e , zc- For this reason, as discussed in |Valiv>] 



ita et al. (2008), all positive interactions with our type of 
interaction are ruled out if Wde is constant. However, the 
situation changes dramatically when we allow for time vary- 
ing Wde- In viable models Wde today is close to —1 and then, 
as —4/5 < w e < 1/3, going towards the past makes Wde less 
negative or even positive. If Wde becomes enough less neg- 
ative before the moment Tde, zc , then the zero crossing can 
actually be avoided. It turns out that the threshold value, 
We,th, depends mildly on all the background density param- 
eters and strongly on the interaction Y / Ho and, of course, 
wo. The larger positive interaction we have, the larger w e we 
need in order to avoid the zero-crossing of pde, and hence 
the singularity of the perturbation equations. This means 
that for a given positive V (and the background parame- 
ters) all the models with w e < w e ,th(Y/Ho, wo) will be miss- 
ing from our Markov chains. There is no similar "top-hat 
cut-off" of models for negative F. This represents a diffi- 
culty in interpreting the marginalized posterior likelihoods. 
Let us assume a completely symmetric situation with re- 
spect to T = 0. Then without the cut-off, we would find 
50% of the area under our Id marginalized posterior for 
Y/Ho to lie in the negative Y region, and 50% in the posi- 
tive T region. However, with the cut-off unavoidably in op- 
eration, even if both positive and negative Y models with 
w e > w e ,th(\Y\/Ho,wo) led to exactly the same theoretical 
predictions (and hence to the same likelihoods), on the posi- 
tive r side the marginalization integral over w e collects only 
the volume 1/3 — w e ,th(Y/ Ho, wo) <C 1.133, while on the 
negative Y side the volume factor is 1/3 — (—4/5) = 1.133. 
The volume factor in the positive Y side becomes smaller 
the larger Y is. Therefore, even in this hypothetical "sym- 
metric" situation, the marginalized likelihood for Y would 
show a strong "preference" for a negative interaction. 



We demonstrate the cut-off effect in Fig. |A1| which 
shows samples from our Markov chains from the run with 
WMAP&BAO data. For example, if Y/H = -0.1, then 
the good-fit region is —0.75 < w e < —0.2, whereas for 
T/Ho = +0.1, all the models -0.8 < w e < -0.45 are for- 
bidden because of the zero crossing of pde- Fig. |A1| shows 
also that when w e is very negative, we can to some extent 
compensate this with less negative wo, as mentioned above. 

With negative values of Y we see in Fig. |A1| another, 
milder, cut-off of models between —4/5 < w e < —2/3. 
Asymptotically with large interactions this cut-off line ap- 
proaches w e = —2/3. From Table [l] on page [3] we can find 
an explanation for this behaviour. Although perturbations 
in the radiation era behave well and we can set adiabatic 
initial conditions, there is a rapidly growing non-adiabatic 
mode in the matter era, if —4/5 < w e < —2/3. This mode 
kicks in faster the stronger the interaction is. In addition, 
this mode grows the faster the further away from w e = —2/3 
we are. Therefore, some models with small interaction rate 
and/or w e close enough to —2/3 survive. 



APPENDIX B: BEST-FITTING MODELS 



We have collected in Table A2 x s and parameters of the 



best-fitting models from our MCMC runs with various data 
sets (Models 1-24), as well as the example models discussed 
in Sec. [2] and Figs. [2] [3J and [9] (Models 25-28). It should 
be noted that Models 6 and 10 (with w e < —0.8) would be 
excluded by the CMB data due to the blow-up of pertur- 
bations, but since the SN and BAO data probe background 
quantities only, these models fit them well. As the x 2 is 
rather insensitive to w e almost equally well-fitting models 
with w e > —0.8 do exist in our chains. 
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Table A2. The x 2 s and parameters of the best-fitting model from each of our 6 MCMC runs for the non-interacting model (T = 0) and 
from our 6 MCMC runs for the interacting model (T / 0). The rows "neg. T" and "pos. T" show the models with largest |T| (found in 
our MCMC chains) such that x 2 is within 4 (~ 2a if the likelihood was Gaussian) from the corresponding best-fitting model. 



Model 


Data 


x 2 


*WMAP 


X 2 <32 


x 2TT 


A/ 2TT 
A^<32 




*BAO 


1 


best r = 


WMAP 


2660.4 


2660.4 


1221.5 


1023.5 


-11.3 


311.4 


22.5 


2 


best T/0 


WMAP 


2659.5 


2659.5 


1221.5 


1022.8 


-11.6 


368.2 


63.7 


3 


neg. r 


WMAP 


2662.6 


2662.6 


1222.0 


1024.9 


-11.1 


389.4 


81.8 


4 


pos. r 


WMAP 


2663.4 


2663.4 


1222.2 


1024.9 


-11.0 


316.1 


27.6 


5 


best r = 


SN 


307.6 










307.6 


1260 


6 


best T/0 


SN 


306.5 










306.5 


146.7 


9 


best r = 


BAO 


3.3 










no a n 


Q Q 
O.O 


10 


best T/0 


BAO 


2.1 










363.7 


2.1 


13 


best r = 


WMAP&SN 


2968.7 


2660.7 


1222.2 


1023.6 


-10.5 


308.2 


7.8 


14 


best T/0 


WMAP&SN 


2968.4 


2660.3 


1221.0 


1024.1 


-11.9 


308.1 


7.8 


15 


neg. T 


WMAP&SN 


2972.1 


2660.8 


1222.2 


1023.6 


-11.9 


311.3 


15.4 


16 


pos. r 


WMAP&SN 


2971.9 


2663.9 


1222.6 


1024.6 


-10.6 


308.0 


13.3 


17 


best r = 


WMAP&BAO 


2667.3 


2660.6 


1222.2 


1023.1 


-10.4 


308.8 


6.9 


18 


best T/0 


WMAP&BAO 


2667.7 


2660.4 


1221.3 


1023.7 


-11.5 


308.1 


7.3 


19 


neg. T 


WMAP&BAO 


2671.7 


2661.4 


1221.7 


1024.2 


-11.7 


309.3 


10.3 


20 


pos. r 


WMAP&BAO 


2671.7 


2664.8 


1224.8 


1023.9 


-7.7 


308.3 


7.0 


21 


best r = 


WMAP&SN&BAO 


2975.7 


2660.6 


1221.8 


1023.5 


-10.7 


308.3 


7.0 


22 


best T/0 


WMAP&SN&BAO 


2975.8 


2660.4 


1220.9 


1024.0 


-11.7 


308.3 


7.1 


23 


neg. T 


WMAP&SN&BAO 


2979.3 


2662.3 


1220.7 


1025.6 


-12.1 


308.7 


8.2 


24 


pos. r 


WMAP&SN&BAO 


2979.6 


2664.3 


1222.7 


1025.5 


-9.9 


308.0 


7.3 




Other models appearing in the figures for phenomenological considerations 








25 


As Model 1, 


but F/H = -0.3 




3680.0 


1223.6 


2026.8 


-9.0 


310.5 


12.3 


26 


As Model 1, 


but T/Hq = +0.3 




3789.9 


1224.5 


2150.0 


-8.8 


312.7 


40.3 


27 


As 25, but uj c (S^deo) & Hq adjusted 


2661.2 


1222.1 


1023.8 


-11.4 


323.3 


41.4 


28 


As 26, but uo c (^deo) & Ho adjusted 


2671.4 


1233.0 


1022.9 


0.2 


307.8 


10.0 



Table A2 - continued The cosmological parameters of the best-fitting models. Today's Hubble parameter Ho is in units kms -1 Mpc -1 , 
the age of the universe is given in Giga years, the distance measure Dy at redshifts z = 0.20 and z = 0.35 as well as the sound horizon 
at decoupling, r s , are in units of h~ x Mpc, where h is defined by Hq = /ikms -1 Mpc -1 . Ampl. denotes the primordial perturbation 
amplitude, indeed ln(10 10 A|). 



Model 






Ho 




r/Ho 


W e 


wo 


n s 


Ampl. 


^deO 


Age 


n 0.20 

V 


n 0.35 


r s 


z dec 


1 


0.0229 


0.107 


61.1 


0.09 





-0.16 


-0.98 


0.975 


3.06 


0.65 


14.2 


552.3 


902.9 


94.2 


1017 


2 


0.0224 


0.133 


49.7 


0.08 


-0.30 


-0.62 


-0.23 


0.954 


3.03 


0.37 


14.6 


507.9 


799.2 


77.2 


1017 


3 


0.0229 


0.196 


47.1 


0.10 


-0.90 


-0.66 


-0.19 


0.971 


3.06 


0.01 


14.7 


500.6 


780.9 


73.2 


1017 


4 


0.0233 


0.091 


58.7 


0.09 


0.12 


-0.13 


-0.81 


0.996 


3.04 


0.67 


14.4 


544.4 


884.4 


91.2 


1018 


5 


0.0757 


0.013 


68.8 


0.08 





0.30 


-0.98 


0.970 


3.04 


0.81 


13.3 


566.2 


937.3 


13.0 


1130 


6 


0.0338 


0.083 


74.6 


0.08 


1.85 


-0.84 


-0.96 


0.970 


3.04 


0.79 


11.3 


565.9 


938.0 


75.5 


1050 


9 


0.0457 


0.010 


68.3 


0.08 





-0.94 


-0.99 


0.970 


3.04 


0.88 


17.4 


582.3 


992.6 


111.7 


1050 


10 


0.0079 


0.010 


89.7 


0.08 


2.92 


-0.87 


-0.98 


0.970 


3.04 


0.98 


12.8 


592.7 


1025 


116.0 


989 


13 


0.0225 


0.109 


67.5 


0.08 





-0.72 


-0.95 


0.960 


3.05 


0.71 


13.9 


559.2 


924.7 


104.0 


1017 


14 


0.0227 


0.116 


67.5 


0.08 


-0.08 


-0.72 


-1.00 


0.963 


3.04 


0.70 


13.9 


560.4 


927.3 


104.3 


1017 


15 


0.0223 


0.133 


62.4 


0.09 


-0.25 


-0.61 


-0.99 


0.955 


3.05 


0.60 


14.1 


550.2 


900.0 


96.8 


1017 


16 


0.0236 


0.084 


64.6 


0.11 


0.19 


-0.06 


-1.00 


1.005 


3.06 


0.74 


14.1 


561.6 


926.3 


99.7 


1019 


17 


0.0226 


0.107 


70.2 


0.09 





-0.78 


-0.99 


0.959 


3.04 


0.74 


13.8 


562.6 


933.3 


106.4 


1017 


18 


0.0226 


0.111 


68.1 


0.08 


-0.04 


-0.66 


-0.99 


0.963 


3.04 


0.71 


13.9 


561.5 


930.0 


105.4 


1017 


19 


0.0220 


0.125 


64.1 


0.09 


-0.21 


-0.70 


-0.97 


0.944 


3.01 


0.64 


14.1 


553.6 


909.8 


100.7 


1015 


20 


0.0228 


0.089 


68.9 


0.09 


0.13 


-0.40 


-0.99 


0.969 


3.03 


0.76 


13.9 


565.1 


938.2 


107.7 


1017 


21 


0.0227 


0.107 


68.8 


0.09 





-0.68 


-0.99 


0.964 


3.05 


0.73 


13.8 


562.6 


933.3 


106.4 


1017 


22 


0.0229 


0.113 


69.1 


0.09 


-0.07 


-0.74 


-1.00 


0.964 


3.05 


0.72 


13.8 


562.6 


933.4 


107.1 


1017 


23 


0.0223 


0.120 


65.9 


0.09 


-0.21 


-0.71 


-0.95 


0.956 


3.01 


0.67 


14.0 


555.5 


915.4 


103.9 


1016 


24 


0.0233 


0.089 


67.9 


0.09 


0.13 


-0.34 


-0.97 


0.983 


3.03 


0.76 


13.9 


562.7 


931.8 


105.8 


1018 




Other models appearing 


in the figures 


for phenomenological considerations 














25 


0.0229 


0.107 


61.1 


0.09 


-0.30 


-0.16 


-0.98 


0.975 


3.06 


0.65 


14.6 


552.9 


905.6 


99.1 


1016 


26 


0.0229 


0.107 


61.1 


0.09 


0.30 


-0.16 


-0.98 


0.975 


3.06 


0.65 


13.9 


551.6 


900.0 


89.1 


1020 


27 


0.0229 


0.137 


56.1 


0.09 


-0.30 


-0.16 


-0.98 


0.975 


3.06 


0.49 


14.4 


538.3 


868.5 


86.6 


1017 


28 


0.0229 


0.080 


66.1 


0.09 


0.30 


-0.16 


-0.98 


0.975 


3.06 


0.76 


14.0 


563.1 


931.1 


102.0 


1017 
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